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Abstract 



It is shown that the well-known non-Abelian static SU(2) black hole solutions have ro- 
tating generalizations, provided that the hypothesis of linearization stability is accepted. 
Surprisingly, this rotating branch has an asymptotically Abelian gauge field with an elec- 
tric charge that cannot vanish, although the non-rotating limit is uncharged. We argue 
that this may be related to our second finding, namely that there are no globally regular 
slowly rotating excitations of the particle-like Bartnik-McKinnon solutions. 



Introduction.— Ever since tlie discovery of a discrete family of particle-like solutions 
of the Einstein- Yang-Mills (EYM) equations by Bartnik and McKinnon as well as 
of their black hole analogues by several workers 0, ||^, it has been a challenge to find 
rotating generalizations of these remarkable objects. This is an exceedingly difficult task. 
It is not even clear how to parametrize the metric, because the coupled field equations 
do not imply the Frobenius integrability conditions for the Killing fields |Q. The usual 
Papapetrou ansatz may therefore be too narrow. Even if this ansatz is adopted, 
the resulting system of partial differential equations becomes considerably more involved 
than in the Abelian case. 

In view of this situation it is reasonable to pursue in a first step a more modest goal, 
based on the assumption of linearization stability Suppose that there is (at least) a 
one-parameter family of stationary black hole or regular solutions of the EYM equations, 
approaching the static solutions mentioned above for angular momentum J = 0, then 
the tangent to this family at J = satisfies the linearized EYM equations. Conversely, 
it is reasonable to expect that for a well-behaved solution of the linearized equations 
around the static configurations there exists an exact one-parameter family of stationary 
solutions. 

In this paper we show that, under the above assumption of linearization stability, 
there are slowly rotating non-Abelian black holes branching off from the (numerically 
known) static ones. To our surprise these black holes are "charged up", in that they 
possess an asymptotically Abelian gauge field with a non-vanishing electric charge. A 
second result, which may also be surprising to many workers in the field, concerns slowly 
rotating globally regular solutions: It turns out that there are no acceptable rotational 
modes in this case; all nontrivial ones develop singularities at the origin. At the end 
of the paper we shall give a heuristic interpretation of how this may be related to the 
charging up of the black holes. We can, of course, not exclude that there is a disconnected 
rapidly rotating branch above some non-vanishing angular momentum. 

Recently, there have been a number of investigations 0, |^ on static axially symmet- 
ric solutions for matter models with YM fields which do not exhibit spherical symmetry, 
but, to our knowledge, we present here the first results on rotational deformations. 
Perturbation equations.— The EYM action for the SU(2) gauge group reads in stan- 
dard notation 

S = I (^--R+-tTF,,F'^^^^d\ (1) 

with Ff^i, = d^Ay — dyA^ + [A^,y4i,]. In the decomposition A^ = A^T^ of the gauge 
potential we choose the gauge group generators T^ = r"/2i with t"' being the Pauli 
matrices. 

An intensively studied family of non-Abelian black hole solutions 0, is described 
by static, spherically symmetric metrics parametrized as 

ds^= (i_'^!!M]a{rfdt^ '^!:L^-r\de' + sm^ edif^), (2) 

\^ r J 1 — 2m[r)/r 

and purely magnetic gauge fields of the form 

A = w{r){-T2 de + Ti sin 9 d(f) + T3 cos 6 dif. (3) 

The functions w{r), m{r), and cr(r) are subject to the EYM field equations. The solutions 
are characterized by the event horizon radius and an integer n counting the number 
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of nodes of w. We shall need the asymptotic behavior of the solutions, 

w = w{rh) + 0{x), m = y + 0(x); w = -1 + ^ + (^-^Y m = M + (^-^^^ , (4) 

at the horizon and infinity, respectively; here x = r — r^, and w{rh), a, and M are fixed 
by the solution parameters Vh and n; the behavior of a is determined by that for w. 
Consider perturbations of a given background equilibrium solution {g^^,A^), 

+ 6gf,^, A^^A^ + 6A^. (5) 

The perturbation equations are obtained by linearization of the field equations: 

SR^, = 2 6T^,, 5(Z}^Fn = 0, (6) 

where D^, = Vy + [A;^, ■ ]. With the notations h^i, = 6g^u and ipfi = ^^fi the perturbation 
equations take the form 

— Vcr V^/i^i. — 2 R^Qu/sh"'^ + W^hau + R^h^n = 4 5T^^, 

- D^D^ij, + Rli,, - 2[F,„, r] + h'^^D^Fp, + V = 0, (7) 

provided that the following gauge conditions are imposed: 

v^K = K = D.r = 0. (8) 

The quantity 5T^y in Eqs. (|^ is obtained by varying the energy-momentum tensor 

T^. = itr {f,^F,p g"^ - \g,.F^^F,,g'^Pg^'^'^ (9) 
with respect to the metric and the gauge field. 

Rotational perturbations.— To identify the most general rotational degrees of free- 
dom, we determine those amplitudes in the partial wave decomposition which can give 
a non-vanishing contribution to the ADM flux integral for the total angular momentum 

= 7^ / e^nk + 5; /l"'^) (fSK (10) 

We do this by expanding h^y and ip^ with respect to an appropriately chosen complex null 
tetrad. (This approach has proved to be extremely efficient for solving wave equations 
0.) The complete separation of the angular variables in Eqs. (|^) is achieved, provided 
that each of the tetrad projections (for ip^ we project in addition onto T± = Ti ± iT2, 
and T3) is chosen as the product of a radial amplitude, depending only on t and r, and a 
spin-weighted spherical harmonic sYim{G^ (p) Here the spin weight is s = 0, ±1, ±2, 
depending on the projection under consideration. Passing back to the coordinate basis, 
we obtain the full mode decomposition for perturbations with respect to the angular 
momentum quantum numbers. (We do not present here the explicit expressions in view 
of their complexity.) 

Since the background solutions are spherically symmetric, group theoretical argu- 
ments imply that modes with different pairs (/,m) decouple, and for a given / (an irre- 
ducible representation of the rotation group) there is a (2/ -|- l)-fold degeneracy labeled 
by m. The angular order m does not occur in the equations for the radial amplitudes. 
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and we can thus choose m = 0. Next, passing to cartesian coordinates x*, we compute 
the ADM angular momentum by integrating over a two-sphere at finite radius r in Eq. 
( [T0| ) and then taking the hmit r — oo. The angular dependence of the integrand implies 
then that the integral vanishes for any r unless / = 1, and that only the ho^ perturbation 
component can give a non-vanishing contribution. Choosing / = 1 and suppressing the 
time-dependence, the tensor modes with the spin weight s = ±2 vanish and Eqs. (|^ 
reduce to a system of 18 coupled equations. 

Now, the transformation behavior of the angular momentum under space and time 
refiections (P,T) implies that only those perturbation amplitudes are relevant which are 
even under P and odd under T. For / = 1 these appear only in Hq^ and in two isotopic 
component of ipQ. They decouple from the remaining modes because the background 
solutions are P and T symmetric. This analysis leads finally to the following most 
general ansatz (up to global coordinate rotations) for the stationary rotational modes: 

h = 2S{r)sm'^edtd(p, i/j = (t^^ sin 9 + cos 6^ dt . (11) 

(The invariance of ip and the background gauge field (^ under P becomes manifest after 
a suitable gauge transformation). Conditions (|^) for the ansatz are fulfilled identically 
and the perturbation equations (|^) reduce to the following coupled system for the radial 
amplitudes S, and r] in ([TT|) : 

-r'Nai — ] + {2N + A '^ ] S + ANr^w' (-)' + ^^"^^ ~ x - v) = 0, 




r^Na I ^ ) +{l+w^- 2w'^N) x-2w7]-rN (w'S)' + \^2Nw'^ + '-j S = 0, 

- w'^n) 7] - Aw x+ '^^^ ~ "^^^ S = 0. (12) 

Here w, m, and a refer to the background solutions, = 1 — 2m/ r. For a solution of 
these equations, the ADM angular momentum is 

f = 5i lim \ . (13) 
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Qualitative considerations.— Before turning to the numerical analysis of the radial 
equations, we make some qualitative remarks. Consider first a Kerr-Newman black 
hole with mass M, angular momentum J and electric charge Q. Let A denote the 
electromagnetic potential of the solution. Next, consider the embedding of this Abelian 
solution into the SU(2) gauge theory, such that the gauge field potential is A = AT^- 
Suppose that \ J\ and \Q\ are small and linearize the solution with respect to J and Q. 
The result can be viewed as a Schwarzschild black hole with two linear hairs. The axial 
gravitational perturbation can be represented in the form of Eq. ( pUf ) with 



S{r) = , 14 

r 

while the Yang-Mills hair is described hy A = T^{Q /r)dt. After a gauge transformation 
with the SU(2)-valued function U = exp((7r — 9)T2) exp{ipT3) this gauge field becomes 

A = ^ UTsU-' + UdU-' = i/j + Ap,,e, (15) 
r 
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where Apure is a pure gauge potential of the form with w{r) = —1, and the pertur- 
bation ip is given by 

^ = ^(Tisin^-T3Cos^)cil (16) 
r 

Now let us return to the perturbation equations (|1^) and concentrate on their 
solutions in the asymptotic region r ^ Vh- In this limit, the background geometry is 
approximately Schwarzschild and w{r) ^ —1, i.e., the background gauge field is almost 
pure gauge. Eqs. (|T^) then split into the two independent groups - one equation for the 
metric perturbation S, and two coupled equations for the Yang-Mills amplitudes x and 
V- 

The acceptable solution of the first equation is given by Eq. ([I^). The equations 
for X and rj admit the solution x(r) = —ri{r) = Q with constant Q, which, together with 
the background pure gauge field, exactly corresponds to the gauge field specified by Eqs. 
([TSp, ([T6|). We therefore conclude that the solution of the perturbation equations in the 



asymptotic region is close to the linearized Kerr-Newman solution with charge Q and 
angular momentum J. 



So far the parameters J and Q are independent. However, Eqs. (^21) split into the 
two independent groups only asymptotically. It is to be expected that the regularity of 
the solution in the entire domain r > Vh will imply relations for the parameters J and Q 
describing the asymptotics. We shall indeed show below that the charge Q of the black 
hole is uniquely fixed by the value of its angular momentum J. 

Numerical solutions.— We now describe briefly the numerical analysis of the full prob- 
lem. Guided by the Abelian example, we are looking for global solutions of Eqs. (|12D 
which are everywhere regular. Note that they do not have to be normalizable, as the 
Kerr-Newman example already shows. Eqs. (p!2D have regular singular points at the 
horizon and at infinity. At the horizon, the formal power-series solution is found to be 
(using Eq. 



S = co + c,x + Oix^), x = -^^^^ + C2X + 0ix'), r] = -^ + C3X + 0{x'), (17) 

rh rh 

where cq, ci, C2, and C3 are four independent integration constants. The solution at 
infinity contains three independent parameters, J, Q, and C4: 

2JM aQ 

J — I :r 

r/ = Q (^-1 + !1_ ) + =^ + O ( ^ ) . (18) 

Clearly, J and Q are, respectively, the angular momentum and the charge of the black 
hole. 

Now, we extend these asymptotic solutions from both sides to the intermediate 
region and impose the matching conditions for the functions Xi a-iid 17, and their 
derivatives. This gives six linear algebraic equations for the seven free parameters in 
([T7|), (rH), and hence the matching can generically be fulfilled. Thus there exists globally 




a solution whose asymptotic behavior is specified by Eqs. (|T^) and (|^). For a given J 
the values of the remaining six coefficients will be proportional to J. In particular, we 
have a relation of the type 

Q = n^{rh) J, (19) 
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Figure 1: The radial behavior of the rotational excitations for the n = 1, rh = 1 
non-Abelian black hole solution. 



where f2„(r/i) is a function determined by the background non-Abelian black hole solu- 
tion. 

The numerical integration of Eqs. ( |12|) (see Fig.l) with the boundary conditions 
specified by Eqs.(p!7|), (0) reveals that in general none of the seven coefficients in {TT}, 
(p!8| ) vanishes. We therefore conclude that, under the assumption of linearization stability, 
each non-Abelian black hole solution admits stationary rotational generalizations with 
electric charge given by Eq. (0). The behavior of the function flnii^h) for the lowest n's 
and for < r/j < cx) is shown in Fig. 2; numerical values are given in Table 1. 



Table 1. Numerical data for n„(r/j). 



n 


rh = 


Th = 0.5 


rh= I 


rh = 5 


1/a = r/, > 1 


1 


-0.4396 


-0.2416 


-0.1117 


-0.0108 


-0.0052 a + 0(^2) 


2 


+0.0092 


-0.0091 


-0.0096 


-0.0003 


-0.0014 a + 0(^2) 


3 


-0.0145 


-0.0084 


-0.0010 


-8 ■ 10-6 


-3- 10-5 a + 0(^2) 



Discussion.— To obtain the asymptotic behavior of n„(r/i) for r/i 3> 1 we make use 
of the fact that, in this limit, the spacetime geometry of the non-Abelian black holes 
is Schwarzschild, up to corrections of order 1/r^. Then, introducing the new radial 
coordinate ^ = ar with a = l/r^, in terms of which the metric functions are N = 
1 — 1/^ + 0(a;^), cr = 1 + O(a^), we expand Eqs. (0) with respect to a. The solution 
in zeroth order of the expansion is S^^^ = —J/^, and rj^'^^ = x^^^ = 0. In first order, we 
obtain two linear equations for 77^^-' and x^^^ with the source terms constructed from S^^\ 
These equations, as Eqs. (|12|), admit solutions which become constant at infinity, thus 
giving rise to an electric charge of order a. 

In the opposite limit, — ^ 0, the functions f2„(r/i) assume finite values. The non- 
Abelian black hole solutions converge for r/j — >■ in the region r > pointwise to the 
regular Bartnik-McKinnon solutions DH, which suggests that the same holds for their 
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Figure 2: The charge-to-angular-momentum ratio Q/J versus the event horizon radius 
for the three lowest branches of rotating black holes solutions. 

perturbations. In other words, one can expect that the particle-like solutions admit 
rotational excitations as well. We have to investigate, however, what happens at the 
origin. There, the behavior of the unperturbed solutions is w = ±(1 — /3r^ -|- 0(r^)), 
m = 2/3^r^ -|- O(r^), f3 being a parameter. Eqs. (|T2|) then give for the perturbations 

5 = -4 ao /3 - ^ (l + 4(3^) r^^ + 02 + 0{r^), 

X = ao (1 - 10/5^2) +a,r+(^a2(3-^y^ + 0{r'), 

ri = ao {1- 8/?2 r^) + ai r + ag + O(r^), (20) 

where oq, ai, a2, and 03 are four integration constants. Notice that the number of free 
parameters, together with those in the asymptotic solution (0), is again large enough 
for matching in the intermediate region. The numerical integration of Eqs. (|T^ for the 
regular backgrounds with the boundary conditions given by ( p^Sf ) and (pO|) shows that 
none of the seven parameters in (|TB|), (^) vanishes, and the solutions in the region r > 
are indeed very close to the black hole rotational modes obtained in the limit ^ 0. 
In particular, the ratio Q/J for the particle-like solutions is given by Q„{0). However, 
although the functions S, x ^^id rj are regular, the corresponding perturbations (pTf ) 
nevertheless become singular at the origin, since the field ansatz ([Tl|) involves the factor 
1/r, while the coefficient ao in Eq. (^) does not vanish. We therefore conclude that the 
regular solutions do not admit rotational excitations. 

The appearance of a nonvanishing electric charge induced by the rotation is a 
somewhat surprising feature of the solutions under consideration - the static non-Abelian 
black holes being neutral. It is natural to wonder where this charge originates from and 
where it resides. These questions cannot be answered in a gauge invariant manner, but 
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Figure 3: The charge function Q(r) given by Eq. (21) for the n = 1, rh ~ 0.22 black hole 



it may be instructive to compute the flux 



7 = i *SF 



(21) 



through a two-sphere of finite radius in the distinguished gauge where the field 5F 
becomes asymptotically Abelian. This gives uniquely 

X 



I V \ r 



I sin^ 9 + 0) cos^ e^dtAdrY (22) 

which motivates us to introduce the charge function 

'r]{r) - 2x{r)\' 



Q(r) = -2tr (T3 /) = 



3r 



(23) 



normalized such that Q(cxo) = Q. The behavior of Q(r), shown in Fig. 3, suggests that 
part of the total charge Q is hidden behind the horizon, and the rest is distributed 
between the horizon and infinity. If approaches zero, the hidden charge assumes a 
finite value, such that its distribution becomes 5-like, making it plausible that there are 
no acceptable rotational perturbations for the solitons. 

To summarize, we have found that there is a branch of slowly rotating colored black 
holes with an asymptotically Abelian gauge field, whose electric flux is proportional to 
the angular momentum J. Mathematically the appearance of the charge is understood 
as follows: The static black holes, similarly to the Schwarzschild solution, admit in the 
far zone two independent hairs, describing rotational and charged excitations. However, 
unlike the situation in the vacuum case, these two hairs become coupled to each other in 
the near zone via the background gauge fleld, and regularity enforces a relation between 
Q and J. Physically, on the other hand, this rotational "charging up" may seem rather 
unusual, since one normally expects that rotation induces only dipole corrections to the 
static field. ~ Our second main result is that the particle-like static Bartnik-McKinnon 
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solutions do not admit continuously connected stationary rotational excitations. Despite 
this fact, we do not see any reason why solitons in other matter models should not have 
rotational states. Rotating magnetic monopoles or Skyrmions might be examples of such 
objects. 
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